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Abstract: Recently, interest has developed in the distribution of interjet energy flows, 
with for example the leading-log calculation of the highly non-trivial colour structure of 
primary emissions in 4-jet systems. Here we point out however that at leading-log level it is 
insufficient to consider only multiple primary emission from the underlying hard antenna — 
additionally, one must take into account the coherent structure of emission from arbitrarily 
complicated ensembles of large-angle soft gluons. Similar considerations apply to certain 
definitions of rapidity gaps based on energy flow. We examine this new class of terms in 
the simpler context of 2-jet events, and discover features that point at novel aspects of the 
QCD dynamics. 
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1. Introduction 

The study of interjet transverse energy (E t ) flow away from jets was suggested by March- 
esini Sz Webber [1] as a method of separating QCD Bremsstrahlung contributions from 
those of the underlying event in hadron-hadron collisions. Recently, in the 2 + 2-jet case 
their considerations on mean energy flows into specified detector regions have been ex- 
tended by Berger, Kucs and Sterman [2] to energy-flow distributions. 

This work apart from its potential phenomenological value at hadron collider exper- 
iments also represents theoretical advances in understanding and computing the complex 
colour topology dictating the flow of soft gluons in 4 jet events, and follows in part from 
earlier work which discussed perturbative calculations for rapidity gaps in terms of interjet 
energy flows [3, 4]. 

One of the main aims of Refs. [2-4] is to resum logarithms L in transverse energy 
at single-logarithmic accuracy i.e. all terms (a s L) n . The authors consider the structure 
of multiple independent soft emissions, and their virtual corrections, from an antenna 
consisting of the four hard partons (the incoming and outgoing jets) and show how it 
exponentiates. This is a rather delicate procedure because of the interferences that arise in 
the colour algebra when squaring the amplitudes involving 4-hard partons and an arbitrary 
number of soft gluons. The above mentioned studies therefore represent a considerable 
advance in the field in terms of understanding the 'colour content' of a multi-parton hard 
scattering (see e.g the discussion in Ref. [3]). 

In certain contexts, for example for the invariant mass distribution of a dijet pair [5], 
this exponentiation of what we shall label as primary emissions is sufficient. However when 
considering energy flows in restricted angular regions precisely as in Refs. [2], or equivalently 
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Figure 1: (a) veto on primary emissions going into Q; (b) veto on energy-ordered secondary 
emission going into f2; (c) veto on emissions going into Q that are radiated coherently from the 
ensemble made up of all (much) harder emissions not going into tt. 

definitions of diffraction based on energy cuts [3,4] (or multiplicity cuts), there is another 
class of single-logarithmic (SL) terms which must be accounted for. To be more specific it 
was shown in Ref. [6] that for any observable that is sensitive to radiation in only a part 
of phase space — a non-global observable — it is necessary to account for secondary and 
higher emissions (as defined below) to obtain an answer correct to SL accuracy. 1 

To illustrate this point better we turn to figure 1, which for simplicity refers to a 2-jet 
event, though the logic remains the same for the case with several hard jets. It shows a 
patch f2 in rapidity and azimuth into which the total (transverse) energy flow is restricted 
to be less than some small amount Qn- There are various sources of SL terms. One source 
comes from vetoing 'primary' emissions which fly directly into f2, figure la, where primary 
means that their radiation pattern corresponds to that for the antenna associated with 
the two hard jets. The lowest-order contribution of this kind comes from an incomplete 
cancellation between real and virtual terms and gives — 2as ^ F An In ^ , where Q is the 
hard scale of the problem and An relates to the area of 0. 

A second source of SL terms comes from diagrams such as figure lb. Here we have 
a large-angle 2 primary emission which flies outside with energy Qi, such that Q S> 
Qi S> Qn — this gives us one power of a s In Forbidding it from radiating a secondary 
emission into Q gives powers of a s ln^-, which after integration over Qi translate into 
a set of SL terms, (a s ln-^) n . This kind of term has been neglected in [2-4], as well as 
in several other contexts [7-9]. To correctly account for it at all orders it is necessary to 
consider soft emission into Q which is coherently radiated from arbitrary ensembles of soft 
(but harder), large-angle energy ordered gluons outside of f2, figure lc, rather than just the 
hard initiating jets in the picture. For simplicity we call this kind of emission a secondary 
emission, though this is only a figure of speech, since it is coherently radiated from external 
ensembles which may consist of primary, secondary, tertiary, etc. gluons. In general we refer 
to the class of terms generated by such contributions as non-global logarithms. 

In earlier work [6] we calculated such terms for observables sensitive to radiation in a 
single hemisphere of a two-jet event. This was done exactly to second order in a s lnJy-, 

1 Another context in which non- global terms would be relevant is if one were to attempt to carry out 
a resummation to single-logarithmic accuracy for certain kinds of isolation criteria for isolated photons. 
However we are not aware of any such resummation being currently in existence. 

2 Strictly speaking, large-angle means of the same order as the angles involved in the definition f2 — if 
Q has a boundary at a small angle to a hard-parton axis, then large-angle can actually mean 'of the same 
order as that small angle.' 
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and numerically in the large- Nq limit at all orders. Here we extend this work to the case 
of interjet energy flow in two-jet events. We view this in part as an intermediate step to 
a calculation in the 3 and 4-jet cases, however it also has value in its own right. Firstly it 
will turn out that an analysis of the dependence of the effect on the geometry of the patch 
VL casts considerable light on the dynamical mechanisms involved in non-global effects. 
Secondly it allows us to make a general order of magnitude estimate of the importance of 
non-global terms relative to those from the resummation of primary emissions. Finally the 
measurement of energy-flow distributions in 2-jet events in e + e~ collisions or DIS could 
well be of intrinsic interest since it would be complementary to measurements in hadron- 
hadron collisions, and in particular, free of the problems associated with the underlying 
event. 



2. Primary emission form factor 

In this paper we shall be considering as our observable the amount of transverse energy E t 
flowing into a patch Q in rapidity and azimuth: 

E t = Y J E t ,i. (2.1) 

We are interested in the probability for E t to be less than some value Qq which is much 
smaller than the hard scale Q of the process in question: 



i rQn ( ] a 

Xn(Qn,Q) = - dE t — , (2.2) 
a J dEt 



where a is the Born-order cross section for the process — in our case the production of 
two jets in e + e~ or of 1 + 1 jets in DIS. 

In order, later on, to quantify the effect of non-global logs it is useful first to calculate 
the contribution to from primary emissions alone. This is the much simpler 2-jet 
analogue of what has been calculated in [2] for 4-jet systems. 

At first order in a s , the logarithmically enhanced contribution to comes from the 
incomplete cancellation of real and virtual contributions for a soft primary emission: 



^\Qn,Q) = -4C F -J -J adn - = -— Aa m — , (2.3) 



a s f Q/2 dk t f j dcf) 4C F a s A Q 
<Qn 

where we have introduced the notation An for area of the region £1, 



= (2.4) 

The upper limit in the kt integral is arbitrary to single-log accuracy, as long as it is of order 
Q. 

When the logarithm of Q/Qn becomes large enough to compensate the smallness of 
a s , it is necessary to include terms (a s In to all orders. If one assumes (incorrectly, as 
we shall see) that multiple wide-angle soft gluons from a two-jet system are simply radiated 
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independently according to a two-particle antenna pattern, then eq. (2.3) can be extended 
to all orders by accounting for the running of the coupling 3 and then exponentiating the 
answer: 

£n,7>(Qn, Q) = (t(Qn, Q)) = exp [-4C F A n t] . (2.5) 

The subscript V on Ef^-p serves as a reminder that we have only taken into account primary 
emissions and t is defined to be the following integral of a s , 

rQ/2 dk t „. _ 1 ,„«s(Q/2) 
'Qn 

where the second equality holds at the one-loop level and (3q = (HCa — 2n/)/(127r). 



t(Qn,Q) = — -rra s (k) = — -In , (2.6) 

27T Jq k t 47r/^o a s (Q n ) 




3. Leading order calculation of non-global effects 

As well as dealing with primary emissions, it is necessary to account also for contributions 
from (secondary) emissions coherently radiated into O from large-angle soft-gluon ensem- 
bles outside of f2. We will denote the contribution from such non-global terms by the 
function S(t), such that to SL accuracy 

Z n (t(Q n ,Q))=S(t)Z niP (t). (3.1) 

To start with, we calculate the leading order contribution to S, i.e. S2, where we define 
the following series expansion for S: 

\ 

\ 

S(t) = ^2s n t n . (3.2) 

n=2 b 

Since this kind of contribution only starts with 
secondary emissions, there is no S\ term. In the 

calculation of S2, we shall be entitled to equate ,*^-Ar|-^\ 
t with gtlngS-. 

The exact value of S 2 depends on the geom- Fi S ure 2: Thc kind of dia S ram to be con - 
r , . , . „ TT . . , sidcrcd for thc calculation of 6>2 in thc case 

etry of the patch \l. Here we calculate it ana- . . ,. ,. „ . , , . 

01 a rapidity slice oi width An. 
lytically for the case where O is a slice in rapidity of width A77. The kind of diagram to 

be considered is shown in figure 2, where a and b are quarks (they may be outgoing or 

incoming depending on whether for example we are dealing with e + e~ or DIS in the Breit 

frame) and 1 and 2 are gluons. We introduce the following four-momenta 

fco = |(l, 0,0,1), (3.3a) 

k b = |(1, 0,0,-1), (3.3b) 

ki = :n^(l,O,sin0i,cos6>i), (3.3c) 

^2 = ^2-^(1, sin 62 sin</>, sin 62 cos (j), cos 9 2 ) , (3.3d) 



3 Strictly speaking the running of the coupling is connected with the collinear branching of the primary 
gluons. This however is a separate issue from that of large-angle soft gluon emission with which we deal 
later on in this article. 
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where we have defined energy fractions x± t2 <C 1 for the two gluons. To our accuracy, we 
can neglect the recoil of the hard particles against the soft gluons. 

We write the squared matrix element for energy-ordered two-gluon emission as (see for 
example [10]) 

w _ Ar (gft) (Ca (al) , C A (bl) ( C A \ (ab) \ 

= C F W X + C F C A W 2 , (3.4b) 

where (ij) = ki ■ kj. The result is valid for 1 3> x\ S> x 2 as well as for the opposite ordering 
of the gluons, and in addition is completely symmetric under interchange of k\ and k 2 . 
(We have however chosen to write it in an asymmetric form so as to emphasise the dipole 
structure of the emissions, namely radiation of gluon k\ from the ab dipole, followed by 
the radiation of gluon k 2 from the al, 16 and ab dipoles). 

The C F piece of the matrix element, W\ corresponds to independent gluon emission 
and is included in the primary emission form factor. To study specifically the modification 
relative to the primary emission case, at this order one must consider the CfCa part of 
the emission probability, W 2 . 

For a general region f2, <S2 is defined through the following equation: 

S 2 In 2 + (ln-§-) = -C F C A [ dcos9 1 < ^-[ dcos0 2 ^ 



->4 />! 



f x 2 dx 2 j x\dx\ 6 [x 2 - W 2 , (3.5) 



{ dcos#i + ^ dcosO^ J dcos6 2 ^ TTT (3-6) 



9L 

16 

which takes into both virtual and real contributions. In the case of a slice of width Ary, 
the angular integrals can be rewritten explicitly 

a a d ^ f a a ^2 
acos0i — — / d cos 8 2 — — 

'2vr 

where we have centred the slice at rj = (the results are independent of its position) and 
defined ±c to be the cosines of the polar angles delimiting the slice, 

Ar ? = lni^. (3.7) 

The integrals over the energy fractions in (3.5) are straightforward. Keeping only the 
leading-logarithmic piece, exploiting the symmetry in Q\ <-> ir — 9\ and performing the 
azimuthal average we have 

/ — C PC 
dcosfli J dcos0 2 F2(cos0i,cos0 2 ) . ( 3 - 8 ) 
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where the angular function F2 is 
F 2 = 



(cos 62 — COS 6\) (1 
Integrating over the polar angles we obtain 



cos #i)(l + cos 62) 



(3.9) 



AC F C A 



12 1 - c 



In — ^— In 




(3.10) 



which can be expressed in terms of A77 as follows: 

g + (Ar?) 2 - At, In (e^» l) - i L i 2 (e^) - i L i 2 (l - , 2 ^) 



S2 = —ACfCa 



where the dilogarithm function is defined as 



U 2 (z) 



-f 

J z 



ln(l - t) 



The functional dependence of £2 on A77 
is shown in figure 3. For small Ar/, S2 goes 
to zero essentially linearly in 77, 

S 2 = -AC F C A [2(1 - ln2Ar/)Ar/ + Ar/ 2 
+0(A V 3 )], (3.13) 

with a logarithmic enhancement due to the 
integrable divergence in eq. (3.8) when 9\ ~ 
62 — —c — thus S2 is roughly proportional 
to the area of the slice. On the other hand 
as Arj increases, S2 rapidly saturates at its 
asymptotic value, 



dt. 



(3.11) 



(3.12) 



o 

u. 

O 
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patch (A0 = At|) 
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-CfCa 



(3.14) 



Figure 3: 1S2 function of 5rj for two dif- 

ferent definitions of ft: a slice in rapidity (using 
eq. (3.11)), and a square patch in rapidity and 
azimuth with 5<p — Sij (S2 determined numeri- 
cally). 



There is a simple physical reason for this 
behaviour: £2 is associated with the dif- 
ference between full coherent emission for 
a pair of gluons, and simple independent 
emission. The dominant contribution to £2 

comes therefore from the region where the two gluons are close together (which by defi- 
nition means the edges of fl since one gluon is in, while the other is out). On the other 
hand, when the two gluons are widely separated in rapidity then independent emission 
becomes a good approximation and there is no contribution to S2 — hence for large Ar], 
1S2 receives no contribution from the centre of the slice, only from its edges, and the value 
of S2 saturates. 
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As well as showing 52 for a slice, figure 3 also shows it (determined numerically) for a 
square patch in rapidity and azimuth. Over the range of Ar] = A(p shown, the behaviour is 
quite different, with an approximate linearity in Ar) — this too can be understood from the 
above arguments: since it is the edges of which contribute dominantly to 1S2, the value 
of 52 will be roughly proportional to the perimeter of the patch, and hence linear in Ar]. 
This holds however only for moderate patch sizes: for very small patches, S2 is roughly 
proportional to the patch area (with a logarithmic enhancement of similar origin to that 
for the slice), while for large patches the periodicity in (ft means that the patch tends to a 
slice. 

4. All-orders treatment 

While for the 2-gluon case the analytical calculations are relatively straightforward, in the 
many-gluon case the situation is vastly more complex: not only does emission into Q come 
from a 'Mexican-cactus'-like external multi-gluon ensemble, but the evolution with t of the 
ensemble outside £1 depends itself on the structure developed at smaller t. Difficulties come 
from both the geometry and the colour structure of the ensemble. As a result, we have not 
succeeded in obtaining even approximate analytical results, and we have to resort to the 
large- Nc approximation and numerical methods in order to extend our calculations to the 
all-orders case. 

4.1 Possible underlying dynamics 

Before going on to the details of the numerical results it is however instructive to consider 
some very rough arguments concerning the underlying dynamics. One way in which Q, can 
stay empty is simply to prevent all emitters outside of f2 from emitting into Q. For each 
emitter 'close' to the edge of ft the price to pay is roughly e~* (we ignore the coefficient of 
t); since, to a first approximation, the typical number of relevant emitters will be of order 
t, we end up with a suppression that goes as e - * 2 . 

Of course when considering very rare configurations, it is risky to base one's argu- 
ments on average properties of the ensemble, such as the typical number of emitters close 
to the edge. For example, instead of suppressing radiation into f2 from gluons in the 
neighbourhood of Q, one could just as well envisage a situation where the neighbour- 
hood of fi is empty, automatically avoiding secondary emissions into £1. We shall try to 
work through this argument with the additional characteristic that we shall discretise the 
problem. Taking Q as a slice of width Ar/ (sufficiently wide that the two edges can be con- 
sidered completely independent), the probability of it staying empty down to some scale 
t is e~ 4CFtAr, S(t). Let us suppose that the condition for secondary radiation not to be 
emitted into f2 is determined by the probability that bands of width Srj ('buffers') on either 
side of the slice stayed empty down to a scale t — St. Then we have a recurrence relation: 

e -AC F tAr, s ^ = e -4C F [tA V +2(t-5t)8 V ] s( t _ ^ j ^ 
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where the term proportional to 5r] in the exponent stems from the suppression of primary 
emission in the 'buffer' bands. Rewriting (4.1) as a differential equation, we have 



5hiS_ 

St 



= -8C F t 



(4.2) 



The factor (Srj/St) can be understood as the average 'speed of motion' (probably propor- 
tional to Ca) of the border between regions with and without emissions {i.e. of the edge 
of the buffer). The resulting form for S is 



S(t)~ e - 4C ^> 



(4.3) 



m 



where we have (arbitrarily) assumed \ jt/ to be independent of t. 

One of the features of this mechanism 
is that resolving only those emissions harder 
than some intermediate scale t\ one will see 
an extended empty buffer region surround- 
ing S7, as in figure 4 (shown for one side of 
a wide slice). The typical size of this buffer 
region should be of the order of 



'buffer 



^buffer - {t 



(4.4) 



Figure 4: A possible structure for emissions 
harder than some intermediate resolution scale 
t' for an event in which f2 stays empty down to 
some scale t > t'. In this scenario, at scale t' 
an empty buffer region exists between f2 and the 
emission closest to £1 



It is the suppression of intermediate-scale 
primary radiation in the buffer region which 
is responsible for the strong suppression of 
S(t) at large t, eq. (4.3). 

An important consequence of the existence of such a buffer region is that the large-t 
behaviour of S(t) would be independent of the shape and size of Q. This is because at 
large t the edges of the buffer region will be far from O and so the details of 0, can have no 
influence on the dynamics at the edge of the buffer. This independence on the geometry of 
Q, together with the explicit observation of a buffer region, would allow us to distinguish 
the buffer mechanism from the alternative mechanism proposed at the beginning of this 
subsection. 



4.2 Numerical results 

We mentioned above that there are two main problems in obtaining all-order results. One 
is the complexity of the colour algebra in the presence of large numbers of gluons. This can 
be eliminated by taking the large- Nq approximation, in which the squared matrix element 
radiation can be broken down into a sum of independent terms each associated with a 
different colour dipole [11]. 

The second source of complexity is the geometry of multi-gluon events, which we 
treat using a Monte Carlo algorithm like that discussed in [6], which essentially models a 
tree-like sequence where a colour dipole emits a gluon, thus branching into two dipoles, 
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each of which may themselves go on to emit (there is a similarity to the Ariadne event 
generator [12]). However unlike a Monte Carlo event generator such as Ariadne, our method 
has the property that it gives results which are purely a function of t (eq. (2.6)), and 
therefore contain just the piece which is guaranteed to be correct, namely the leading 
logarithmic piece. 

The Monte Carlo algorithm returns a function <Smc(*) m the large- Nc limit. We choose 
to correct this function so that at least at order aig the result is correct beyond the large- Nc 
approximation. Accordingly, in what follows, we shall consider 

2Cp 

S(t) = [S MC (t)} Ca , (4.5) 



rather than <Smc itself. 




0.1 0.2 0.3 0.4 0.5 0.6 0.1 0.2 0.3 0.4 0.5 0.6 

t t 



Figure 5: A comparison between (a) the full result for S and (b) the simple exponentiation of 1S2. 
Beyond t = 0.5 the differences between the curves for the finite and infinite slices are not significant 
compared to the statistical and systematic errors involved in the determination of S. 

In figure 5a we show the function S(t) for three different geometries for Q: a square 
patch in rapidity and azimuth (A77 = Acfi = 2), a finite slice in rapidity (A77 = 1) and 
an infinite slice in rapidity. We extend the t scale beyond the phenomenologically relevant 
region in order to better illustrate the general features of S in the different cases. The right- 
hand plot, figure 5b shows what would be obtained if there were a simple exponentiation 
of the S2 term, S ~ exp(<S2i 2 ). 

There are various points to be noted: firstly, as t increases, S decreases with a behaviour 
roughly consistent with a Gaussian suppression, at least for t < 0.5. Secondly, at large t, 
though the normalisation of S(t) depends on the geometry of Q, its behaviour in t seems 
to be universal. This is to be compared to the geometry dependence that would be present 
if the all-orders result stemmed from a simple exponentiation of the £2 term, fig. 5b. We 
note that up to t ~ 0.5, there is a strong similarity between the actual t-dependence and 
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that observed in the exponentiation of S2 for an infinite slice. Beyond t ~ 0.5 however the 
suppression in the full calculation seems to be grow much faster. 

0.25<t<0.35 0.55<t<0.65 




0.1 0.2 0.3 0.1 0.2 0.3 0.4 0.5 0.6 



Figure 6: Contour plots for the distribution of the size of the buffer region, ^buffer, as a function 
of the intermediate resolution scale V ' . For the larger t range, the irregularities of the contours are 
an artifact due to limited statistics. The darkest contour corresponds to --j-^ 2 — = 2. 

The fact that the t-dependence of S is the same regardless of the geometry of seems 
to suggest that the 'buffer' mechanism postulated above might well be responsible for the 
all-orders behaviour of >S. To examine this idea in more detail one can use the Monte 
Carlo simulation to establish whether there really is a buffer region, and to examine its 
size, ^buffer; as a function of an intermediate resolution scale t' < t. Figure 6 shows contour 
plots for the distribution of //buffer values as a function of t' . Here Q is an infinite slice, and 
t is defined as the scale of the hardest secondary emission in O originating from coherent 
emission off partons which are to the right of Q. Correspondingly the buffer region being 
considered is that to the right of f2. 

For the lower t range the buffer region tends to be fairly small — for example at t' = 0.1, 
the most likely buffer size is r/buflfer — 0.7. However, increasing t to 0.55 < t < 0.65, the 
most likely size of the buffer region at t' = 0.1 grows to ?7buffcr — 3.5. Such an increase 
in buffer size for fixed t' as one increases t is precisely what one would expect from the 
'buffer' mechanism postulated earlier, eq. (4.4). From that equation one can also evaluate 
an effective value for (5rj/5t) ~ 9. Using eq. (4.3), one then finds that the coefficient of t 2 
in the exponent is roughly twice 1S2 for an infinite slice. 

This is not quantitatively consistent with what is seen in figure 5, however several points 
should be borne in mind: (a) the arguments in section 4.1 are based on average properties, 
whereas in the end we are interested in the properties of rare events. So for example 
there is a certain spread in the distribution of ^buffer an d this can affect the quantitative 
predictions by a non-trivial factor, as can the fact that 6rj/6t is itself only defined in an 
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average sense (b) one's estimate for {Srj/St) depends somewhat on how exactly one deduces 
it from the plots, e.g. whether as a derivative with respect to if or t; this is connected to 
point (c), namely that our assumption of a constant {8rj/St) (independent of t') may be 
an oversimplification. Indeed if (Sr]/5t) is constant then eq. (4.4) one would expect the 
centre of the distribution of r/buffcr to be depend linearly on t' . However in fig. 6 there 
seems to be some non- linear dependence of the typical ^buffer on t', though it is not clear 
whether this is an artifact of t not being sufficiently asymptotic, or whether there is extra 
dynamics which remains to be taken into account, such as a pile-up of emissions close to 
the edge of the buffer. Furthermore at t ~ 0.5 the behaviour of S{t) departs from the initial 
approximate Gaussian, and starts to fall much faster, also suggesting either that t < 0.5 is 
not yet asymptotic, or that (Sr]/St) has some extra t' dependence. 

In conclusion, we believe that while we may have understood some of the gross features 
of the dynamics involved in the all-orders behaviour of S(t), in particular the importance of 
a buffer mechanism in determining the geometry-independent large-t suppression of S(t), 
there are some significant details which remain to be understood. 
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Figure 7: Xn(t) for two different definitions of O: in the left-hand figure we consider the energy 
flowing into a slice of width Ai] = 1, while in the right-hand figure we consider the energy flowing 
into a square patch in rapidity and azimuth, of size At? = A</> = 2. The upper plots show the 
relation between t and Qn for different values of the centre-of-mass energy Q. 



To understand the phenomenological significance of the non-global logarithms, it is 
interesting to compare results for T,(t(Qn,Q)) with only primary emissions and with the 
full non- global treatment. This is done in figure 7 for two different geometries of Q, a slice 
in rapidity and a square patch in rapidity and azimuth. 
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The upper plots show how t is related to Qa, for different values of the centre of 
mass energy Q. Even at a very energetic next-linear-collider, optimistically trusting the 
calculation down to Qq = 0.5 GeV one only just goes beyond t = 0.25. For current energies 
the largest value of t that is accessible is in the range of 0.15 to 0.2. Taking t = 0.15 as 
our reference value, the inclusion of the non-global effects increases the suppression (on a 
logarithmic scale) relative to that for just primary emissions by a factor of between 1.5 
(patch) and 1.65 (slice). To state it a different way, ignoring non-global effects at t = 0.15 
overestimates the cross section by between 30% (patch) and 65% (slice). At larger t values, 
these figures rapidly become even more dramatic. 

So the non-global effects are not only important from the point of view of the formal 
correctness of the leading-log series, but also numerically significant. 

As an aside we comment on the feature of a divergence in the distribution of Qn at the 
Landau pole that has been observed in certain circumstances in [2,3]. Such a behaviour 
arises when the suppression from En,p(t) is not sufficient to compensate the divergence 
in dt/dQn, i.e. when < 1. However the inclusion of the non-global factor S(t) 

ensures that £n(t) always goes to zero much faster than dt/dQn diverges. Therefore the 
distribution of Qq should go to zero at the Landau pole regardless of the size of C ^ Q " ■ 

Finally, we note that for a comparison to data it would also be necessary to take 
into account non-perturbative effects. One way of doing so might be in terms of power 
corrections [13]. For the mean Qn, in a normalisation in which the power correction to the 
e + e~ thrust has a coefficient c T = 2, the coefficient for a region Q would be cq = AqQ 
where An is the area of Q. For the corresponding differential distribution one should 
roughly expect a shift of the distribution by the same amount. 

5. Conclusions 

As has already been pointed out in [6], in the resummation of any observable sensitive 
only to emissions in a limited angular region of phase space, there is a class of single- 
logs — 'non-global logs' — which leads to a breakdown of the picture of independent 
primary emissions and strict angular ordering. That picture has been quite widely adopted 
in the literature, in certain instances wrongly [2-4, 7-9], at least to the accuracy that 
was claimed. 4 To deal with this class of terms it is necessary to consider emission from 
arbitrarily complicated ensembles of energy-ordered large-angle gluons lying outside the 
region of sensitivity of one's observable. 

In this article we have considered non-global logarithms for observables such as the 
distribution of energy flow in restricted angular regions between jets. This kind of mea- 
surement was originally advocated, in [1,2], for 2 + 2-jet events in hadron-hadron collisions. 
Here we have studied a simpler case, that of 2- jet events, in order to concentrate on the 
specific features of these non-global logarithms, without the difficulties that arise from the 

4 There exist other cases, e.g. [14,15], whose final results are for global variables, but where intermediate 
steps of the derivation make reference to non-global hemisphere- variables without the inclusion of non-global 
logs. It is important not to use those results outside the context of the specific derivation for which they 
are intended without taking care of non-global effects as required. 
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complex colour structure of four-jet systems. We also point out that experimental studies 
of interjet energy flows in e + e _ and DIS, insofar as they are free from contamination by 
the underlying event, may provide useful complementary information to that which could 
be extracted in hadron-hadron colliders. 

Our studies have led us to various conclusions. Most importantly perhaps, is that non- 
global effects, as well as being formally of the same order as those from primary emissions, 
are also numerically almost as large. 

More academically, one of the interesting aspects of non-global logs that is specific 
to interjet energy flows, as opposed to the single-hemisphere observables studied in [6], is 
that they depend on the geometry of 0, the region in which one measures the energy flow. 
At second order in a s (the first order at which these effects appear), for moderately sized 
patches, the magnitude of non-global effects is roughly proportional to the perimeter of 
the patch. In contrast, effects due to primary emissions scale as the area of the patch. 

At all orders, we have only a large-Ac* numerical treatment for non-global effects. 
Despite this it is possible to obtain some interesting insights into the mechanisms that are 
of relevance. One notable result is that, modulo a geometry-dependent normalisation, the 
asymptotic t-dependence of the non-global suppression factor S seems to be independent of 
the size and shape of f2. This fact, together with figure 6, lends support to the hypothesis 
that at intermediate scales there is an empty 'buffer region' around O and that a significant 
part of the non-global suppression factor S stems from the suppression of intermediate- 
scale primary radiation in the buffer region. This allows us to postulate that in other 
processes, such as 2 + 2-jet production at hadron colliders, one will see similar results, with 
the asymptotic t-dependence of S being given by that of figure 5a, raised to a power which 
depends on the number of jets and whether they are quark or gluon jets. 
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